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One of the most serious problems of the staggered fermion lattice action is flavor symmetry violation. Smeared 
link staggered fermions can improve flavor symmetry by an order of magnitude relative to the standard thin link 
action. Over the last few years different smearing transformations have been proposed, both with perturbatively 
and non-perturbatively determined coefficients. What hindered the acceptance and more general use of smeared 
link fermions until now is the relative difficulty of dynamical simulations and the lack of perturbative calculations 
with these actions. In both areas there have been significant improvement lately, that I will review in this paper. 


1. INTRODUCTION 

Fermionic actions with smeared link gauge con¬ 
nections have gained popularity in recent years. 
For staggered fermions, where the effect of smear¬ 
ing on flavor symmetry restoration is easy to un¬ 
derstand perturbatively and straightforward to 
measure numerically, smearing became almost 
mandatory |^,||,^|^. For Wilson-like fermions, 
including overlap and domain wall fermions as 
well, the acceptance of smearing has been slower, 
though there is evidence that smearing improves 
the chiral symmetry of Wilson fermions and in¬ 
creases computational efficiency for all formula¬ 
tions The lack of interest in smeared 

link Wilson fermions is partially due to the fact 
that there is no well defined physical quantity, like 
flavor symmetry violation for staggered fermions, 
that would indicate the effect of smearing for Wil¬ 
son like fermions. 

Smearing can have undesirable effects as it can 
distort the (lattice) short distance properties of 
the configurations |^. What level this distortion 
is depends on the specific smearing and is clearly 
observable form the short distance part of the 
heavy quark potential. While excessive smear¬ 
ing is never desirable, for certain operators, like 
heavy fermions, only very little or no smearing is 
acceptable. Fortunately there is no need to use 
the same smearing transformation, or even the 
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same fermionic formulations, for heavy and light 
quarks, so this problem is easily solvable. 

Smearing both for staggered and Wilson 
fermions can be considered as part of a systematic 
improvement program where irrelevant operators 
are added to the action to remove lattice arti¬ 
facts. Perturbative improvement systematically 
removes lattice corrections of order 0 (a") and 
0 ( 5 "), while the fixed point action program re¬ 
moves lattice artifacts using a non-perturbative, 
renormalization group based approach. At the 
end both methods modify the fermionic action 
by including non-nearest neighbor quadratic and 
higher order terms and modify the gauge connec¬ 
tion by replacing the simple link by some sort 
of extended paths. The systematic fixed point 
action program has been used mainly with Wil¬ 
son like fermions while the systematic pertur¬ 
bative improvement program has been applied 
for staggered fermions, but both methods can 
(and to some extent have been) extended to both 
fermionic formulation. 

In this talk I will review some of the recent 
developments regarding staggered fermions. The 
two smeared staggered actions that are used in 
dynamical simulations at present are the so called 
Asqtad and HYP actions. Since the Asqtad ac¬ 
tion has been reviewed in detail last year |l|, I 
will concentrate on the newer HYP action [|10[ . 
Since the HYP action contains SU(3) projected 
gauge links, it cannot be simulated with the usual 
molecular dynamics based algorithms. I will re- 
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view an updating method that is based on the 
well-known pseudo-fermion update and is efficient 
for smeared link fermions . Finally I show 

some preliminary results obtained with this algo¬ 
rithm for HYP staggered fermions. 

2. THE ADVANTAGES OF SMEARED 
LINK ACTIONS 

I start with a brief outline of the perturba¬ 
tive and non-perturbative methods of construct¬ 
ing smeared link staggered fermion actions. The 
perturbative method is systematic, in principle it 
can be extended to arbitrary order in 0(a") and 
0{g'^), though at present only a fully 0{a?) im¬ 
proved action, the Asqtad action, has been used 
in numerical simulations ||r^ , |l^p^ . The non- 
perturbative method I will review is the hyper- 
cubic smearing (HYP) transformation |^. At 
present it is not part of a systematic improvement 
program, yet it is a smearing transformation that 
appears to be about a factor of two better in fla¬ 
vor symmetry restoration than the Asqtad action 
and incidentally it is also 0(a?) improved. 

Perturbative smeared link actions: Stag¬ 
gered fermions break flavor symmetry as gluons 
at the edge of the Brillouin zone connect fermions 
with different flavors. The tree level pertur¬ 
bative improvement program removes the flavor 
changing gluons by requiring that the smeared 
gauge link vanishes at momenta (tt, 0,0,0), 
(tt, 7r,0, 0), etc. This can be achieved by the ad¬ 
dition of three different kind of gauge loops of 
length three, five, and seven, connecting nearest 
neighbor fermions 



This smearing is frequently referred to as Fat7 
fl^ . The construction cannot remove all flavor 
changing gluons, those with momentum near but 
not exactly tt do not vanish, but their coefficients 
are reduced. Removing the smeared gluons at the 
corners of the Brillouin zone introduces an O(a^) 
term in the small momentum gluon. This term 
can be removed by the addition of yet another 


gauge term, the so called Lepage term j^. To 
make the action fully 0{a^) improved the pure 
gauge part of the action is chosen to be the tad¬ 
pole improved 1-loop Symanzik action and the 
Dirac operator is modified by the addition of a 
third neighbor Naik term, constructed with the 
original thin links, leading to the Asqtad action. 

At tree level perturbation theory SU(3) pro¬ 
jected (unitarized) links cannot be distinguished 
from the non-projected ones. By choice, in the 
Asqtad action the extended gauge loop connec¬ 
tions appear linearly, no projection to SU(3) is 
included. The Asqtad action has been studied 
extensively by the MILC collaboration and many 
of its properties have been reviewed last year in 



Recently a more ambitious program of 1-loop 
improvement has been developed and preliminary 
results were reported at this conference. The 
1-loop, Oia^g^) improved action contains about 
a dozen new four fermion coupling terms, some 
with imaginary coefficients. In numerical sim¬ 
ulations the four fermion terms can be resolved 
by the introduction of auxiliary scalar fields, but 
the presence of imaginary couplings could create 
more problems. To date the action has not been 
used neither in quenched nor in dynamical simu¬ 
lations. Most of the work so far concentrated on 
identifying the different terms of the action, de¬ 
termining their coefficients and estimating their 
importance. It would make future numerical sim¬ 
ulations easier if only a few of the theoretically 
necessary terms were important, especially if the 
imaginary terms could be neglected. The anal¬ 
ysis of the perturbative action therefore is very 
important. 

Non-perturbative smeared link actions: 

To design a non-perturbative improvement 
scheme first one has to understand the origin 
of flavor symmetry violation non-perturbatively. 
Staggered fermions break flavor symmetry be¬ 
cause each Dirac and flavor component of the 
fermion field occupies a different lattice site of the 
hypercube, and couples to different gauge fields. 
If the gauge links within a hypercube fluctuate 
strongly, the different fermionic components ob¬ 
serve very different gauge environment and act 
like independent fields, and not like components 










3 


of a Dirac spinor, thus causing flavor symmetry 
violation. This very simple physical picture sug¬ 
gests that in order to improve the flavor sym¬ 
metry of staggered fermions one has to minimize 
the gauge held fluctuations at the hypercubic or 
plaquette level. More precisely, one should min¬ 
imize the tail of the plaquette distribution, i.e. 
maximize the smallest plaquettes of the configu¬ 
rations, as these are the ones causing most of the 
flavor symmetry violations. In a quenched study 
with different smearing transformations the re¬ 
lation between the minimum plaquette distribu¬ 
tion and flavor symmetry violation was clearly 
demonstrated suggesting a fast, non-perturbative 
method to optimize smearing [T^ . 

HYP smearing was constructed to be very local 
and optimized through the minimum plaquette 
distribution. It is a succession of three levels of 
modified SU(3) projected APE smearing steps, 
each level modified such that at the end none of 
the original links extend beyond the hypercubes 
that attach to the smeared link 


= Projsu{z)[{^-oii)Ui,i^ 
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In Eq. d denotes the HYP smeared links, 
are the original thin links, while and 
are intermediate smeared links. There are 
three parameters in HYP smearing, the APE 
parameters of each level. They can be optimized 
non-perturbatively by minimizing fluctuations at 
the plaquette level. One can also optimize the pa¬ 
rameters of the HYP links perturbatively by re¬ 
quiring that the smeared gauge links vanish at 
the edges of the Brillouin zone. The perturbative 
condition leads to parameters that are, at least 
after tadpole improvement, consistent with the 


non-perturbatively optimized parameter values. 
In numerical simulations the non-perturbatively 
optimized parameters 

ai = 0.75, a2 = 0.60, as = 0.35 (2) 

are used. 

HYP smeared links are fairly local. That is best 
seen in the measurement of the static quark po¬ 
tential. The static potentials measured with thin 
gauge links and with HYP smeared links agree, 
up to an additive constant, at distances r/a > 2 
IM • The difference between the potential values 
is significant only at r/a = 1, and even there it 
can be well described by the perturbative lattice 
Coulomb potential. An added bonus of the HYP 
potential is that its statistical errors are about an 
order of magnitude smaller than the thin link po¬ 
tential errors, making the measurement of the po¬ 
tential parameters, especially the string tension, 
feasible even on relatively small data sets. 

The tail of the plaquette distribution that is 
maximized in HYP smearing contains mainly 
those plaquettes, frequently called dislocations, 
that are responsible for the residual chiral symme¬ 
try breaking in domain wall fermions, and for the 
large real eigenmodes that slow down the over¬ 
lap fermion computations. A smeared link ac¬ 
tion that is optimized to remove these modes will 
likely improve domain wall and overlap computa¬ 
tions as well Q. 


2.1. Flavor Symmetry Restoration in 
Smeared Link Actions 

The effect of smearing on flavor symmetry is 
illustrated in Figure where the relative split be¬ 
tween the Goldstone pion mo and lightest non- 
Goldstone pion ttItt 


A. 


- mg 

mo 


( 3 ) 


is plotted as the function of the square of the 
Goldstone pion mass measured in units of rg 
fr^ . The triangles, circles and filled squares in 
Figure are results obtained with the standard 
thin link, Asqtad and HYP fermion actions on 
8^ X 24, /3 = 5.7 Wilson plaquette gauge configu¬ 
rations with lattice spacing a ~ 0.17fm. The open 
squares correspond to dynamical HYP configura¬ 
tions at similar lattice spacing. They agree with 
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Figure 1. Flavor symmetry violation with 
the standard thin (triangles), Asqtad (circles) 
and HYP (filled squares) staggered fermions on 
a = 0.17fm quenched configurations. The open 
squares correspond to dynamical HYP data at 
similar lattice spacing. The vertical line at 
(romc)^ = 1 correspond to me ~ 400MeV pi- 
ons. 


the quenched HYP data, indicating that improve¬ 
ment seen in quenched simulations carry over to 
dynamical ones. For all mass values the HYP ac¬ 
tion has about an order of magnitude smaller fla¬ 
vor symmetry violations than the standard thin 
link action, and it is about a factor of two better 
than Asqtad fermions. It is an interesting ques¬ 
tion why HYP fermions have better flavor symme¬ 
try restoration than Asqtad fermions. Both are 
perturbative tree level improved, but the HYP 
links have three levels of SU(3) projections and 
more gauge loops at each level. Whether it is the 
SU(3) projection, the extra gauge loops, or both, 
that make HYP fermions better is not obvious. 
Recently it has been shown jl^ that actions with 
projected SU(3) links that have the same tree 
level 0{a^) improvement are also equivalent at 
1-loop level. If that perturbative statement holds 
true non-perturbatively as well, one could argue 
that the SU(3) projection is the key to the im¬ 
proved flavor symmetry of the HYP action. The 
non-perturbative assumption can be easily tested 
either by quenched measurements of flavor sym¬ 
metry with projected Asqtad/Fat7 actions or by 


simply measuring the distribution of the tail of 
the plaquette with these actions. 


2.2. Perturbative Properties of Smeared 
Actions 

A frequent criticism toward smeared (or, for 
that matter, any kind of improved) actions is that 
their perturbative properties are not known. This 
excuse has always been a weak one, and now, 
that several papers have been published show¬ 
ing that perturbative calculations with smeared 
links, even projected ones, are doable, the criti¬ 
cism now longer holds. Rather, one should ask 
if the perturbative properties of smeared actions 
are better or worse than the standard thin link ac¬ 
tions. One can recall that with thin link actions, 
both with staggered and Wilson fermions, the 1- 
loop perturbative matching factors for most phys¬ 
ically interesting quantities are large. That makes 
the connection between the lattice and continuum 
schemes like MS unreliable. It has been observed 
that in the case of staggered fermions the pertur¬ 
bative contribution from the fermion doublers is 
as large as from the usual gluon tadpoles [l^ . 
Since smearing removes the fermion doublers, it 
should also improve the perturbative matching 
factors. Several recent works showed that this is 
indeed what happens [15 20,1^,^. For example 


in the 1-loop matching coefficients for several 
staggered bilinear operators were calculated for 
different smeared actions. These matching coeffi¬ 
cients connect continuum to lattice as 


^cont ^ + (4) 

-l-ofT — diagonal terms. 


If a coefficient Ci^i « 5 at 1/a = 2GeV lattice 
spacing, the 1-loop perturbative corrections are 
about 10%. Table Q shows the scalar-scalar coef¬ 
ficient that is relevant for the fermionic conden¬ 
sate for the standard thin link, Asqtad and HYP 
actions. Obviously both the Asqtad and HYP ac¬ 
tions satisfy the minimum requirement Ci^ < 5, 
their matching factors are 10-30 times smaller 
than the thin link value. The HYP value is even 
small enough to use on rougher, 1/a = 1.2GeV 












5 


action 

Thin link 

Asqtad 

HYP 

Cl,l 

-29.3 

- 2.2 

- 0.6 


Table 1 

The perturbative matching factors for the scalar- 
scalar bilinear for the standard thin link, Asqtad 
and HYP fermions pO| , pT[ 

configurations where one finds 

E(2GeH)^ = ^Eiatt(l/a = 1.2GeV), (5) 

^p/s 

with 

Zp,s = 0.92. ( 6 ) 

This perturbatively determined renormalization 
Z factor agrees with the value obtained via the 
non-perturbative matching method of I will 
discuss that further in Sect. 4. 

3. DYNAMICAL SIMULATIONS WITH 
SMEARED ACTIONS 

Smeared actions that are linear in the thin 
gauge links, like the Asqtad action, can be sim¬ 
ulated with the usual molecular dynamics based 
algorithms. Calculating the fermionic force for 
these actions, while tedious, is straightforward. 
The only problem observed with these simula¬ 
tions is the unusually large autocorrelation time 
for the topological charge . 

Actions with projected smeared links, like the 
HYP action, can not be simulated using molec¬ 
ular dynamics methods as the fermionic force 
is prohibitively expensive to evaluate. However 
these actions are smooth enough in the gauge 
fields to be simulated using a modified version 
of the well known pseudo-fermion algorithm, the 
partial-global stochastic metropolis or PGSM up¬ 
date ||^,^, |rT|jl^ . The PGSM algorithm satis¬ 
fies the detailed balance condition [|4|, and with 
smeared link fermions it can be an efficient algo¬ 
rithm. 

We consider a smeared link action of the form 

S=Sg(U) + Sg{y)+Sf{V), (7) 


where Sg{U) and Sg{V) are pure gauge actions 
depending on the thin links {[/} and smeared 
links {U}, respectively, and S'/ is the fermionic 
action depending on the smeared links only.The 
smeared links are constructed deterministically, 
like hypercubic (HYP) blocking. Any gauge ac¬ 
tion can be used for Sg{U) while Sg{V) is cho¬ 
sen to optimize the algorithm. The fermionic ac¬ 
tion describing n / degenerate flavors of staggered 
fermions is 

Sf{V) = -^\nnetn{V) ( 8 ) 

with M defined on even sites only. In 

the PGSM update first a subset of the thin links 
{U} are updated and a new thin gauge link con¬ 
figuration {U'} is proposed with transition prob¬ 
ability that satisfies detailed balance with Sg(U). 
The proposed configuration is accepted with the 
probability 

Pace = min{l,e-^®«e-«*[^-^l«}, (9) 

where A5g = SgiV ) — Sg{V) is the dif¬ 
ference in the smeared gauge actions, A = 

H “"/Z®, and the stochastic vector 

^ is generated with Gaussian distribution. In Eq. 
^only one stochastic vector is used in each update 
step, the acceptance probability approximates the 
ratio of the old and new actions 

Det > (10) 

only on average. If the stochastic formula esti¬ 
mates this ratio poorly, the autocorrelation time 
of the update will be large, it could even be infi¬ 
nite. This can be seen if we consider the standard 
deviation of the stochastic estimator 

det “^(A) =< >^, ( 11 ) 

ai = < >2 

= det “^( 2 A — 1 ) — det ( 12 ) 

Eq. ^is valid only if the matrix 2A —1 is positive 
definite. If the matrix A has even one eigenvalue 
that is less than or equal to 1 / 2 , the formula in 
Eq. 1^ is not valid, the standard deviation is 
infinite. Unfortunately this is frequently the case 
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as the eigenvalues of the fermionic matrices and 
n vary from l/4m^ to above 16. 

There are several ways to improve the efficiency 
of the stochastic estimator |Q. Here I mention 
only one, the determinant breakup. If we rewrite 
the fermionic action of Eq. as 

=-In Det (13) 

with rib an arbitrary integer, the fermionic deter¬ 
minant can be estimated using nbnfjA random 
source vectors as 

det =< 

with A = 

standard deviation is finite as long as none of the 
eigenvalues of A is less than 1/2"'’, a condition 
that is much easier to satisfy. In fact for every 
finite quark mass it is possible to choose rib such 
that the condition is satisfied. The determinant 
breakup reduces the statistical fluctuations of the 
stochastic estimator by taking the sum of f /4 
smaller terms instead of the original term 

ni,nf/4 

-r[H-i]c^- X] (15) 

i=l 

While taking the average of several stochastic 
estimates in the acceptance step of the orig¬ 
inal PGSM algorithm (Eq. ^) would violate 
the detailed balance condition, the determinant 
breakup procedure is still exact. 

An added bonus of the determinant breakup is 
that now simulating arbitrary number of degen¬ 
erate or non-degenerate flavors is straightforward 
as long as rib is chosen to be a multiple of 4. 

With the determinant breakup of Eq. the 
acceptance probability can be close to the the¬ 
oretical maximum, . The 

PGSM update will be effective only if this value 
is close to one, i.e. the proposed conhgurations 
{!/} —> {[/'} in the first step of the algorithm are 
close to the dynamical configurations of S. Fairly 
good agreement can be achieved if we choose 
Sg{U) such that it matches the average plaquette 
and/or correlation length of the dynamical con¬ 
figurations. That requires the shift of the pure 
gauge coupling that can be compensated by the 


second pure gauge action term, Sg{V). Since this 
term depends on the smeared links, its fluctua¬ 
tions are greatly suppressed and including this 
term in the acceptance rate does not change it 
signihcantly. 

The fractional powers of the fermionic matrices 
in A can be evaluated using polynomial approx¬ 
imation but to evaluate the acceptance 

probability according to Eq. requires nbUf/A 
multiplication with A, a considerable computa¬ 
tional overhead. A detailed study of the auto¬ 
correlation of the PGSM update can be found 
in iQ. Figure || shows the cost, measured in 
matrix multiplies, of creating statistically 
independent configurations that are separated by 
two autocorrelation time lengths as the function 
of the number of links changed in {U} to cre¬ 
ate {U } in the first step of the PGSM algo¬ 
rithm. The data was obtained with rif = 2 flavors 
of dynamical HYP fermions on 8^ x 24, about 
lOfm^ volume configurations with different de¬ 
terminant breakup parameters {rib = 8,12) and 
quark masses corresponding to mT^/rrip « 0.7 and 
0.55. In case of the heavier quark mass the opti¬ 
mal update touches about 4000 links or around 
10% of the total links on this volume. With 
an acceptance rate of 50%, 2 x 10® M mul¬ 
tiplies are needed to create statistically indepen¬ 
dent configurations. For comparison, creating in¬ 
dependent configurations with the standard thin 
link staggered action at similar lattice spacing at 
TO,r/wp ~ 0.7 would require about 6 — 8 x 10® 
matrix multiplies, a few times more than with 
the HYP action using the PGSM algorithm. 

The PGSM algorithm is a volume square al¬ 
gorithm. Greating independent configurations on 
larger volumes will become more and more expen¬ 
sive. However, judging from Figure the com¬ 
putational cost of creating even 100-200fm^ con¬ 
figurations is not more than 5-10 times more ex¬ 
pensive than molecular dynamics update of thin 
link staggered fermions. 

4. SIMULATION RESULTS WITH HYP 
FERMIONS 

The results summarized in this section are pre¬ 
liminary, obtained on modest size lattices at large 
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Figure 2. Cost of creating independent configu¬ 
rations with n/ = 2 flavor HYP action, measured 
in fermionic matrix multiplies. Results are from 
lOfm"*^ lattices with lattice spacing a ~ 0.17fm. 
Open squares: rib = 8 , mT^/mp « 0.7, filled 
square: rib = 12, m-„/mp « 0.7, crosses: rib = 12, 
rriT^Imp « 0.55 simulations. 


lattice spacing with limited computational re¬ 
sources. They are in no way comparable to the 
extensive data collected with the Asqtad action 
over the last couple of years, some of it reported 
at this conference |^,^. My goal in this section 
is to show that the PGSM method is applicable, 
and at least for standard quantities, nothing un¬ 
expected is seen with the HYP action. 

Simulations with the HYP action were per¬ 
formed on 8^ X 24 lattices with a « 0.17fm lattice 
spacing with three quark masses corresponding to 
roughly rriT^/mp = 0.7, 0.65 and 0.55. To moni¬ 
tor finite size effects and to study the scaling of 
the algorithm the simulations were repeated on 
12^ X 36 lattices with the two smaller quark mass 
values. While no autocorrelation time measure¬ 
ment is available on the larger lattices, other fac¬ 
tors indicate that the algorithm scales with the 
square of the volume, maybe even better when 
comparing the small 8^ x 24 lattices to the larger 
volume. 

In Figure I have already showed that flavor 
symmetry restoration with the dynamical HYP 
action is consistent with the quenched result. Fig¬ 
ure H shows the scaling of the Goldstone pion 



Figure 3. The Goldstone pion mass squared as 
the function of the quark mass, in units of rg. 
Triangles: 8 ^ x 24, squares: 12^ x 36 lattices, all 
with rg ^ 3.0. 

mass squared with the quark mass, both mea¬ 
sured in units of rg. The large finite volume dis¬ 
tortion on the smaller volumes is surprising at 
first, quenched simulations with similar rriT^Imp 
ratios and volumes do not show such effect. How¬ 
ever the data is actually consistent with finite vol¬ 
ume chiral perturbative calculations p8|j29| . Ac¬ 
cording to these calculations finite volume cor¬ 
rections of the data above are expected to be 
a few percent on the smaller volume with the 
largest quark mass and also on the larger volume 
with the smallest quark mass, larger for the other 
data points on the smaller volume. Using data 
with small finite volume corrections we can esti¬ 
mate the renormalization factor Zpjg based on 
the non-perturbative method of [^. The value 
we obtain is consistent with the perturbative pre¬ 
diction in Sect. 2.2 

■^p/s = 0-92(1) (non — perturbative). (16) 

Gan we see chiral logarithms in this data? Fig¬ 
ure U shows the ratio / rriq as the function of 
the quark mass, again in units of rg. The data 
points with small finite volume corrections are 
consistent with a constant value, at these quark 
masses there is no obvious sign for chiral loga¬ 
rithms. We can use this data to predict the chi¬ 
ral condensate E through the GMOR relation. 
With the renormalization factor obtained above 









Figure 4. m^/m as the function of the quark 
mass, in units of tq. Notation is the same as in 
Figure ||. 

the chiral condensate, translated to MS scheme 
at 2GeV is predicted to be 

^{2GeV)rl = 0.36(4). (17) 

This value, though with large errors, seems 
to be larger than the quenched value 
S^^"(2Gey)rQ = 0.30(4) predicted by several 
groups using a much more sophisticated finite 
volume method that can be applied only with 
exactly chiral fermions Hi!®- The difference 
between the quenched value and Eq. 0 could 
be only numerical accuracy, though one would 
expect a larger condensate from dynamical sim¬ 
ulations. New results at smaller lattice spacing 
and larger volumes can resolve this question in 
the near future. 

I conclude this section with some recent re¬ 
sults on the topological susceptibility with HYP 
fermions. At small quark masses the topological 
susceptibility is expected to scale linearly with the 
square of the Goldstone pion mass with a coeffi¬ 
cient that is predictable from chiral perturbation 
theory (solid line in Figure]^). Higher order cor¬ 
rections might push the value even lower (dashed 
curve in Figure I) [I- Numerical simulations 
have a hard time reproducing this behavior, es¬ 
pecially at small quark masses and large lattice 
spacing, even with improved fermions |^,^. In 
Figure ^ I show results that were obtained with 
thin link staggered fermions at a ~ O.lfm lattice 


Figure 5. The topological susceptibility as the 
function of the Goldstone pion mass squared, in 
units of tq. Triangles: standard thin link stag¬ 
gered fermions at o « 0.1 Ofm |Q; open squares: 
HYP action on 8^ x 24 lattices, a ~ 0.17fm; 
filled squares: HYP action on 12^ x 24 lattices, 
a w 0.17fm. The dashed curve is the prediction 
from , the short lines on the right indicate the 
quenched value. 


spacing and with HYP fermions at a « 0.17fm 
lattice spacing. The topological charge was mea¬ 
sured after 2-3 levels of HYP smearing using an 
improved topological charge operator Q. The 
difference between the smaller and larger volume 
results with HYP fermions is mainly due to the 
change of rg and m,r due to finite volume effects. 
The topological susceptibility is quite indepen¬ 
dent of the volume. While even HYP fermions do 
not reproduce the theoretical expectations, they 
seem to be consistent with the thin link action re¬ 
sults at 70% smaller lattice spacing. It is unlikely 
that the difference between the numerical results 
and theoretical expectations are due to the rem¬ 
nant flavor symmetry violations. Even the heav¬ 
iest of the non-Goldstone pions would predict a 
smaller topological susceptibility than the mea¬ 
sured value at the smallest quark mass. The reso¬ 
lution to this problem might come from exploring 
the difference of what objects the fermions iden¬ 
tify as instantons versus what objects the gluonic 
charge operators identify as instantons in these 
simulations. 
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5. CONCLUSION 

Smeared gauge links remove short scale vac¬ 
uum fluctuations, improving the flavor symme¬ 
try of staggered fermion actions. In addition 
to improved flavor symmetry, smeared fermion 
actions have improved perturbative properties, 
like 10-30 times reduced 1-loop matching fac¬ 
tors. While these actions are more complicated 
than thin link ones, with the newly developed 
partial-global stochastic Metropolis (PGSM) up¬ 
date even SU(3) projected actions can be simu¬ 
lated on fairly large volumes. There are already 
extensive numerical data with the dynamical Asq- 
tad action using the molecular dynamics based R 
algorithm. The HYP action has been used in ex¬ 
ploratory studies only with the PGSM algorithm. 

Many of the properties of the smeared stag¬ 
gered actions should apply to Wilson like fermions 
as well. Smearing could improve chiral symmetry, 
and reduce the occurrence of exceptional config¬ 
urations. Smeared actions also require smaller 
clover coefficients. In domain wall fermions 
smearing can reduce the residual chiral symme¬ 
try breaking while for overlap fermions smearing 
increases computational efficiency. These proper¬ 
ties are worth exploring in the future. 
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